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1.  Introduction 

The  orbits  of  over  500  comets  have  been  calculated  from  observational  data 
obtained  while  the  comets  were  visible  from  the  earth.  From  these  observations 
it  is  found  that  the  orbits  fall  into  two  main  groups: 

(1)  Those  with  parabolic  or  nearly  parabolic  orbits. 

(2)  Those  that  have  definitely  elliptic  orbits. 

The  second  group  contains  all  comets  with  periods  less  than  100  years.  It  is 
found  that  of  over  70  comets  that  go  to  make  up  this  group  only  three  have  an 
inclination  to  the  ecliptic  of  more  than  30°  and  the  average  inclination  is  13°. 
Also,  apart  from  Halley’s  comet,  all  have  the  same  direction  of  motion  around 
the  sun  as  Jupiter.  These  two  facts  lead  to  the  conjecture  that  originally  all  these 
comets  belonged  to  the  first  group  and  have,  by  one  or  several  close  encounters 
with  Jupiter,  been  diverted  into  their  present  orbits.  Comet  Brooks  (1889  V) 
was  observed  to  undergo  this  capture  by  Jupiter. 

Both  Lyttleton’s  [1]  and  Oort’s  [2]  theories  on  the  origin  of  the  comets  assume 
that  the  short-period  comets  were  forced  into  their  present  orbits  by  the  pertur¬ 
bative  effects  of  the  planets  (in  particular  Jupiter)  and  that  originally  these 
comets  were  in  osculating  parabolic  long-period  orbits. 

Van  Woerkom  [3]  investigated  the  perturbations  of  comets  which  came  near 
to  Jupiter  and  he  found  that  the  root  mean  square  values  of  the  change  in  the 
reciprocal  of  the  semimajor  axis  of  comets  with  perihelion  distances  of  1  and 
4.5  a.u.  was  78  and  128  a.u.-1  X  10-6  respectively.  The  problem  attempted  in 
this  paper  is  an  extension  of  Van  Woerkom’s  work.  By  a  suitable  choice  of 
parameters  the  effect  of  Jupiter  is  evaluated  for  all  parabolic  comets  whose  dis¬ 
tance  of  closest  approach  to  the  sun  lies  in  the  range  0.04  to  4  a.u.  The  lower  limit 
of  0.04  a.u.  is  taken  because  Lyttleton’s  accretion  theory  postulates  a  stream  of 
sun-grazing  comets. 

Jupiter’s  effect  can  be  measured  in  terms  of  the  change  in  the  reciprocal  of 
the  semimajor  axis  of  the  orbit  and  an  estimate  of  the  distribution  of  this  quantity 
is  made  for  each  value  of  the  perihelion  that  is  used. 
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2.  Calculation  of  the  perturbations  due  to  Jupiter 

All  the  calculations  are  based  on  the  energy  changes  as  the  comet  describes 
a  parabolic  orbit. 

If  E  is  the  energy  per  unit  mass  of  a  comet  in  describing  a  parabolic  orbit  about 
the  sun,  then  the  change  in  the  energy  8E  for  a  complete  orbit  will  be  calculated. 
A  sufficient  gain  in  energy  to  make  E  positive  will  mean  that  the  comet  would 
escape  from  the  solar  system,  and  a  negative  E  means  that  it  will  return  in  an 
elliptic  orbit.  In  calculating  the  change  in  energy  we  will  assume  that  Jupiter 
is  the  only  planet  that  has  a  significant  effect. 


m 


(x.*v) 


The  sun  and  Jupiter  are  assumed  to  be  rotating  about  their  common  center 
of  gravity  G  with  constant  angular  velocity  n.  At  any  time  t  the  comet  is  at  a 
distance  r  from  the  sun  and  at  a  distance  p  from  Jupiter  (ST  is  the  direction  of 
Jupiter  at  t  =  0).  The  energy  E  of  the  comet  can  then  be  obtained  from  the 
following  considerations  (using  figure  1) : 

The  comet  m  is  at  ( x ,  y,  z)  (rotating  axes  as  shown),  so  its  equations  of  motion 


are 
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(2.1) 

x  —  2  ny  —  n?x  =  G  7-  | 
J  dx 

(f+0 

(2.2) 

Q 

ij  +  2  nx  —  n2y  —  G  —  \ 
dy 

(f+9 

and 

(2.3) 

It 

<r  p) 

where  G  is  the  gravitational  constant  and  S  and  J  are  the  masses  of  the  sun  and 
Jupiter  respectively. 

If  we  multiply  (2.1)  by  x,  (2.2)  by  y ,  (2.3)  by  z,  integrate,  and  then  add  to¬ 
gether,  we  obtain  Jacobi's  energy  integral 

(2.4)  |  (»  +  +  i’-)  -  \  +  f)  =  G  (^  +  j)  +  C, 

where  C  is  a  constant. 

If,  instead,  m  is  referred  to  fixed  axes  at  G,  with  X  in  the  direction  of  T  and 
Z  coinciding  with  z,  then  the  above  integral  becomes 

(2.5)  |  (P  +  P  +  Z2)  -  n(XY  -  XY )  =  G  +  C, 

since  X  =  x  —  ny  and  Y  =  y  +  nx.  The  last  integral  may  be  written, 

(2.6)  E  —  nh3  —  constant, 

where  E  —  (l/2)u2  —  G(S/r  +  J/p)  =  total  energy,  and  /i»  =  XY  —  XY  = 
angular  momentum  about  GZ. 

Therefore,  to  find  8E  (the  change  in  energy),  it  will  be  sufficient  to  find  8hz. 
The  couple  at  G  has  its  third  component, 

C  J 

(2.7)  ht  =  y  (XY j  -  YXj)  +  ^  (XYS  -  YXS), 

where  the  coordinates  of  the  sun  are  (Xs,  Y s,  0)  and  those  of  Jupiter  are 

(Xj,Yj,0). 

Now  let  r,  a,  /3  be  the  spherical  polar  coordinates  of  m  relative  to  parallel 
axes  at  S,  that  is,  pole  in  the  direction  SZ,  (3  measured  from  S  T;  then 

(2.8)  X  +  ds  cos  nt  =  r  sin  a  cos  /3, 

(2.9)  Y  -j-  ds  sin  nt  =  r  sin  a  sin  /3, 

(2.10)  Z  —  r  cos  a. 

T  he  coordinates  of  J  and  S  are 


(2.11) 


Xj  =  dj  cos  nt, 
Xs  =  —  ds  cos  nt, 


Yj  =  dj  sin  nt, 

Ys  =  —  ds  sin  nt, 


Zj  =  0, 
Zs  =  0. 
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Therefore, 

(2.12)  XY  j  —  YXj  =  rdj  sin  a  sin  ( nt  —  p) 
and 

(2.13)  XY s  —  YXS  =  —rds  sin  a  sin  (nt  —  /3); 
and  since  Sds  =  Jdj  —  n,  say,  the  required  component  hi  is 


(2.14)  hi  =  nG  y—  —  r  sin  a  sin  (nt  —  fi). 

Z 


WHEN  m  IS  AT  At 

Figure  2 

Comet’s  path  about  the  sun. 

Now  suppose  m  is  describing  a  parabola  about  S,  which  is  not  in  the  plane  of 
SJ  motion,  then  Ai  is  the  perihelion;  AXN  is  the  great  circle  through  Ax  which 
meets  the  ecliptic  in  N  the  ascending  node,  and 

T N  =  longitude  of  node  =  9, 

7r  —  TNAi  =  inclination  =  i, 

(2.15) 

NAi  =  angular  distance  of  perihelion  from  node  =  w, 

AiPi  =  v. 

Note  the  direction  ST  is  the  direction  of  Jupiter  at  time  t  =  0.  The  direction 
cosines  of  SP\  give 
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cos  8  cos  (a>  +  v)  —  sin  8  sin  (a;  +  t>)  cos  i  =  sin  a  cos  /3, 

(2.16)  sin  8  cos  (a>  +  v)  +  cos  6  sin  (a>  +  v)  cos  i  =  sin  a  sin  0, 

sin  (w  +  y)  sin  t  =  cos  a. 

Hence, 

(2.17)  sin  a:  cos  (nt  —  13) 

=  cos  (a)  +  v )  sin  (nt  —  9)  —  sin  (a>  +  y)  cos  i  cos  (nt  —  0) ; 

therefore 

(2.18)  8E  =  7z5/i3 

=  nG»  ^ 

where 

(2.19) 

/(0  =  cos  (co  4-  v )  sin  (nt  —  6)  —  sin  (co  +  v)  cos  (nt  —  6)  cos  i, 
p2  =  d2  +  r2  —  2c?r[cos  (co  +  v)  cos  (nt  —  6)  +  sin  (a>  +  v )  sin  (nt  —  8)  cos  i], 

and  d  =  ds  dj  (that  is,  d  is  the  distance  between  the  sun  and  Jupiter). 

The  equation  of  the  parabola  that  m  is  describing  about  S  is 

(2.20)  r  =  q  sec2 1> 


where  q  =  SAi  (that  is,  distance  of  closest  approach). 

It  should  be  noted  that  if  the  comet  were  describing  an  ellipse  or  a  hyperbola, 
then  the  equation  of  the  orbit  would  contain  an  extra  parameter,  that  is,  the 
eccentricity  of  the  orbit. 

If  we  let  r  =  tan  v/2,  then  (2.20)  becomes 


(2.21) 


r  =  q(  1  +  T2), 

1  v 

dr  —  -  sec2  -  dv. 


Now,  for  the  undisturbed  parabolic  orbit  about  the  sun  we  have  by  Kepler’s 
laws 

(2.22)  <  =  (f)'7(T+^ 

=(177 (i+r,)dT 
=(!7  VH 

if  we  take  r  =  0  when  t  =  0. 

Also,  since  Jupiter  is  in  a  circular  orbit  about  the  sun,  we  have  S/d 2  =  n2dj. 
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Thus  (S  -f  J)/d  —  n2d2,  since  S/dj  =  J/ds  =  (S  +  J)/d.  Therefore,  n2  = 
C 8  +  J)/d\ 

Let 


(2.23) 


_  (2*s  +  A1/2; 

\s  d*  ) 

but  since  J  is  negligible  in  comparison  with  S  we  may  make  the  approximation 


(2.24)  X  =  (f)m 

It  is  normal  to  express  the  perturbative  effects  not  in  terms  of  the  energy 
change  but  in  terms  of  the  change  in  the  reciprocal  of  the  semimajor  axis  of 
the  comet’s  orbit,  that  is,  8(1/ a).  Since  the  binding  energy  of  a  nearly  parabolic 
orbit  is  GS/2a,  where  G  is  the  gravitational  constant,  S  is  the  mass  of  the  sun, 
and  a  is  the  semimajor  axis  of  the  comet’s  orbit,  we  have 


(2.25) 


where 

(2.26) 


SJ 

M  -  (iS  +  j)d 


Sds  —  Jdj. 


3.  Numerical  evaluation  of  5(1  /a) 

To  evaluate  8(1 /a)  we  break  the  integral  into  two  parts, 


(3.1) 

and 

(3.2) 


_  2 nn  C‘ 

"  s  . 


m 


dt 


Ij  = 


— 2wjj 

S 


/: 


rM. 


dt. 


If  we  make  the  substitution  nt  =  X[r  +  (l/3)r*]  in  Is  and  use  the  relation¬ 
ships  r  =  <?(1  +  r2)  and  r  =  tan  v/2,  we  get 

/„-#?/■ 


(3.3)  Is  -  j*,  (1  +  dr’ 

where 

(3.4)  /i(r)  =  (1  —  r2)  cos  nt( cos  w  sin  6  +  sin  w  cos  6  cos  i) 

+  2 t  sin  nt( sin  ca  cos  6  +  cos  w  sin  0  cos  t). 


PERTURBATIONS  OF  COMETARY  ORBITS  155 


Hence  the  parameters  a>,  0,  and  i  can  be  separated  out  and  the  problem  is  now 
reduced  to  evaluating  the  following  two  integrals, 

(3.5)  I\  =  Imag  jf  ^  j  exp  (t  +  |  r3^  J  dr, 
and 

(3.6)  h  =  Real  J  exp  |\*X  +  |  ^  J  dr. 

If  the  path  of  integration  is  taken  as  the  contour  r  =  r  exp  (71-/64)  the  integrands 
die  away  very  rapidly  and  the  evaluation  is  trivial.  These  integrals  were  eval¬ 
uated  on  the  Manchester  University  Ferranti  Mercury  Computer  for  the  fol¬ 
lowing  values  of  X: 

(3.7)  X  =  0.001  (0.001)  0.01  (0.01)  0.1  (0.1)  1 
(when  X  =  0.001,  q  '=.  0.04  a.u.). 

The  method  of  evaluation  used  was  that  of  Gaussian  quadrature,  which  is 
incorporated  in  a  standard  library  program  for  infinite  integrals  [4]. 

The  integration  of  Ij,  however,  is  not  as  simple  as  this  because  of  the  presence 
of  branch  points  of  the  integrand  in  the  complex  r  plane.  If  r  =  u  +  iv  then  it 
can  be  shown  that  for  large  r  these  branch  points  lie  on  the  curves 


(3.8) 


v  = 


1 

±X(1  +  u2) 


log 


and  their  positions  on  these  curves  are  given  by 


(3.9) 


u  =. 


3(2 n  ±  1)tt  +  4> 
X 


where  </>  and  k  are  constants  which  depend  on  the  parameters  of  the  comet’s  orbit. 

Because  of  these  branch  points  it  is  impossible  to  deform  the  contour  as 
was  done  in  the  case  of  Is- 
If  we  make  the  substitutions 


nt  —  x  = 


(3.10) 


then 

(3.11) 

where  /2(x)  =  /i(r) 
and 


t  =  tan|> 

r  =  q(  1  +  r2), 
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(3.12)  g2(x)  =  (b  1  +  6(1  -f  r2)2  —  2  (cos  ( x  —  0)[(1  —  r2)  cos  <o  —  2 r  sin  co] 

+  cos  i  sin  ( x  —  0)[(1  —  r2)  sin  co  +  2r  cos  co]})3/2, 

where  b  =  X/V2. 

It  can  be  shown  that  for  sufficiently  large  x0 


(3.13) 


9i(x) 


(s)'7. 


sin  co  cos  i  cos  x  —  cos  co  sin  x 


/\\b!z  r  *  gjn  w  gjn  x  _j_  cog  w  cog  x  cog  i 

+  \v  l - ** - rfx 


(3.14)  f°^dx  =  0(x  or). 

J  X  0 

Therefore  the  remaining  terms  of  the  asymptotic  expansion  are  O(xo2). 

It  will  also  be  noticed  that  the  sufficiently  large  x0  the  terms  in  the  asymptotic 
expansion  behave  like 


<*■“>  />&' 

The  method  that  suggests  itself  from  this  is  that  we  integrate  directly  as  far 
out  as  a-o  and  then  from  xQ  onward  we  integrate  in  steps  of  t  and  the  contribu¬ 
tions  from  each  step  are  added  to  the  value  of  the  integral  so  far  obtained.  In 
this  way  an  alternating  sequence  is  formed  and  the  limit  of  the  sequence  can 
be  estimated  by  the  standard  Euler  method.  It  was  found  that  for  values  of 
a-0  >  20tt,  the  integrand  showed  sufficient  oscillatory  character  for  the  Euler 
technique  to  be  used. 

If  8(1 /a)  is  considered  as  a  function  of  the  parameters  6  and  co  then  the  fol¬ 
lowing  relationships  hold, 

5  (a)  ^  =  8  (a)  [0  +  T>  w  +  *L 

(3.16) 

8  (a)  [tt  -  9,  r  -  co], 

and  when  i  =  0  and  tt,  the  parameters  9  and  co  can  be  combined  as  one  parameter 
<f>.  Therefore  if  a  step  of  r/8  is  taken  in  each  of  the  parameters  i,  6,  and  co,  the 
following  values  are  sufficient  to  cover  the  sphere:  i  =  'k/8(tt/8)7tt/8,  6  = 
0(ir/8)4ir/8,  co  =  0(7r/8)157r/8;  and  when  i  =  0  and  -it,  <f>  =  0(x/8)7-ir/8.  The 
values  of  X  used  are  as  in  the  evaluation  of  Is,  where  X  =  (2<?3/d3)1/2,  and  where 
q  is  the  perihelion  distance  and  d  is  the  distance  between  the  sun  and  Jupiter. 
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0 


Figure  3 

Graph  of  |5(l/o)p|  as  a  function  of 
for  the  case  perihelion  distance,  q  =  4.130  a.u., 
that  is,  for  Jupiter  X  =  1. 

4.  Results 

In  calculating  the  energy  change  it  was  assumed  that  the  effect  of  the  other 
planets  was  small  when  compared  with  Jupiter’s  effect.  To  show  this,  similar 
calculations  to  the  above  were  performed  for  Saturn,  Uranus,  and  Neptune  and 
the  energy  changes  produced  by  these  four  major  planets  were  compared.  The 


Graph  of  Az  as  a  function  of  4> 
for  the  case  i  =  0  and  X  =  0.01,  0.1  and  10, 


158 


FOURTH  BERKELEY  SYMPOSIUM:  KERR 


Graph  of  Az  as  a  function  of  w 
for  the  case  X  =  0.001  and  i  =  0,  ir/4  and  ir/2. 

Upper  panel :  0  =  0;  middle  panel:  9  =  x/4;  lower  panel:  9  =  *■/ 2. 
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q  ina.u. 

Figure  6 

Graph  of  Az  as  a  function  of  q 
for  the  case  «  =  0,  9  =  ir/4.  and  ir/2. 

Upper  panel :  i  =  a-/ 4;  lower  panel :  i  =  a-/2. 

quantity  5(l/a)p  was  evaluated  for  theu  for  major  planets  for  i  =  0,  4>  =  0(ir/8)ir, 
and  for  a  perihelion  distance  of  4.130  a.u.  (this  corresponds  to  X  =  1  in  the  case 
of  Jupiter).  The  i  =  0  plane  was  chosen  because  there  is  always  a  singularity 
of  Ij  in  this  plane  and  hence  8(1 /a)p  will  always  obtain  its  maximum,  that  is, 
8(1 /a)p  will  become  infinite  when  the  comet  collides  with  the  planet. 
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The  results  are  set  out  in  graphical  form  in  figure  3  and  for  convenience 
\8(l/a)p\  was  plotted  against  <£.  It  will  be  seen  that  the  effect  of  Uranus  and 
Neptune  is  negligible  except  when  the  comet  is  nearly  in  collision  with  these 
planets.  However,  Saturn's  effect,  while  smaller  than  Jupiter’s,  is  significant,  and 
in  ignoring  it  we  have  introduced  errors  which  can  be  large  in  the  case  of  comets 
which  pass  close  to  Saturn.  In  calculating  these  curves  it  was  assumed  that  the 
the  planets  moved  in  circular  orbits. 

Typical  results  are  set  out  in  figures  4,  5,  and  6.  For  convenience  the  substitu¬ 
tion  A z  =  10B5(l/a)  has  been  made.  From  these  it  will  be  seen  that  not  only 
does  \Az\  depend  on  how  close  the  comet  comes  to  Jupiter  but  also  on  the  angle 
at  which  the  planes  cut,  that  is,  as  i  tends  to  7r/2  the  magnitude  A z  decreases. 
This  remark,  however,  is  not  always  true  because  there  are  singularities  of  Ij 
other  than  in  the  i  —  0  plane  (for  example,  there  is  one  at  approximately 
i  —  ir/2,  6  =  7ir/8,  co  =  10ir/8,  and  X  =  1),  and  the  effect  of  these  will  tend  to 
invalidate  any  general  comments. 

6.  Calculation  on  particular  comets 

The  program  was  also  used  to  evaluate  A z  for  six  of  the  comets  that  have  been 
observed  to  be  hyperbolic  at  their  perihelion.  The  effects  of  both  Jupiter  and 
Saturn  were  taken  into  account.  In  making  these  calculations  it  was  assumed 
that  the  comets  were  in  parabolic  orbits,  and  that  the  planets  were  moving  in 
circular  orbits  around  the  sun.  Galibina  [5]  has  also  calculated  A z  for  20  comets 
and  where  the  two  lists  overlap  her  results  have  been  included  in  table  I  (which 
gives  the  post  orbit)  for  comparison.  From  these  results  it  will  be  noted  that 
only  Comet  1914  V  will  return. 


TABLE  I 


Comet 

Az 

units  10~6  a.u.-1 

Az 

(Galibina) 

1914  V 

3.69 

3.8 

1922  II 

-58.86 

— 

1925  I 

-63.12 

-59.1 

1925  VII 

-32.55 

— 

1932  VI 

-27.56 

-26.1 

1936  I 

-35.95 

-49.1 

Van  Biesbroeck  [6]  has  also  calculated  As  for  Comet  1914  V  and  he  obtained 
a  value  of  12.6. 


6.  Distribution  of  Az[Az  =  1056(l/a)] 

If  it  is  supposed  that  the  parabolic  comets  for  a  given  X  (that  is,  a  given  peri¬ 
helion  distance  q )  have  their  orbits  disposed  at  random  then  we  require  the 


PERTURBATIONS  OF  COMETARY  ORBITS 


161 


distribution  function  for  Az.  Consider  the  coordinates  (i,  d)  which  determine  the 
plane  of  the  comet’s  orbit.  Then  if  dS  is  an  element  of  area  on  the  celestial 
sphere  we  have 

(6.1)  dS  =  C  sin  i  dd  di, 

where  C  is  a  normalizing  factor  which  depends  on  the  range  of  variables. 
Therefore,  since  co  is  uniformly  distributed,  the  distribution  of  A z  is  given  by 

(6.2)  C  sin  i  di  dd  du  =  —C  d{ cos  i)  dd  da>. 

This  suggests  that  we  ought  to  make  cos  i,  d,  and  w  independently  uniformly 
distributed.  Therefore  in  the  above  calculations  it  is  necessary  to  weight  the 
values  of  A z  associated  wTith  each  i  in  such  a  way  that  the  distribution  of  cos  i 
is  uniform. 

Now  A z  has  been  evaluated  for  nine  values  of  i  given  by  ir  where  ir  =  rir/8 
with  r  =  0(1)8;  then  the  weight  to  be  associated  with  each  of  these  values  of  ir  is 

TABLE  II 

Values  of  q  and  <t  Corresponding 
to  28  Values  of  X 


X 

q  (in  a.u.) 

<r 

0.001 

0.0413 

77.9 

0.002 

0.0656 

77.2 

0.003 

0.0859 

77.8 

0.004 

0.1041 

78.7 

0.005 

0.1207 

80.7 

0.000 

0.1364 

81.2 

0.007 

0.1511 

81.3 

0.008 

0.1652 

79.2 

0.009 

0.1787 

78.2 

0.010 

0.1917 

77.7 

0.020 

0.3043 

79.2 

0.030 

0.3987 

75.8 

0.040 

0.4830 

75.7 

0.050 

0.5605 

76.5 

0.060 

0.6329 

78.9 

0.070 

0.7014 

806 

0.080 

0.7667 

77.3 

0.090 

0.8293 

75.0 

0.100 

0.8897 

73.1 

0.200 

1.4123 

68.9 

0.300 

1.8506 

60.6 

0.400 

2.2418 

55.5 

0.500 

2.6014 

52.6 

0.600 

2.9376 

47.7 

0.700 

3.2556 

46.8 

0.800 

3.5587 

50.1 

0.900 

3.8494 

47.8 

1.000 

4.1295 

49.3 
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(6.3) 


COS  IfX]  COS  — i 
2 


r  5^0  or  8, 


COS  i\  —  COS  2-o 


2 


r  =  0  or  8. 


Therefore  each  triple  A z(i,  6,  u>)  is  given  a  weight  according  to  the  value  of 
i  associated  with  it.  The  r.m.s.  value  A z  is  therefore  the  square  root  of 
]£u>(Az)2/]£w,  where  =  128. 

Now  let  this  r.m.s.  value  be  a;  then  table  II  gives  the  values  of  <r  and  q  cor¬ 
responding  to  the  28  values  of  X.  These  values  are  also  plotted  against  q  in  figure 
7  and  a  straight  line,  in  the  sense  of  least  squares,  is  fitted.  The  line  obtained  is 

(6.4)  a  =  —9.24  q  +  80.8 

and  the  sum  of  the  squares  of  the  residuals  =  263. 

If  this  formula  is  applied  to  calculate  a  for  the  two  cases  given  by  Van 
Woerkom  [3],  that  is,  q  —  4.5  a.u.  and  q  =  1  a.u.,  we  obtain  the  following  results, 

a  —  42  when  q  =  4.5  a.u., 

(6.5) 

a  =  72  when  q  =  1  a.u., 


which  compare  with  Van  Woerkom ’s  results  of 

a  =  128  when  q  =  4.5  a.u., 

(6.6) 

a  =  78  when  q  =  1  a.u. 

The  distribution  of  A z  is  found  by  summing  the  weights  of  all  the  triples 
A z(i,  6,  co)  which  give  values  of  A z  in  each  of  the  intervals  —300(10)300  in  units 
of  10_5(a.u.)_1.  The  contribution  from  each  interval  is  then  divided  by  128  to  give 


Figure  7 

Plot  of  <r  as  a  function  of  q. 

The  fitted  line  is  a  =  — 9.24g  +  80.8. 
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Figure  8 

Comparison  of  double-exponential, 
Gaussian,  and  empirical  distribution  of  Az. 
Upper  panel:  q  =  0.0413  a.u.,  a  =  77.9; 
middle  panel:  q  =  1.8506  a.u.,  a  —  60.6; 
lower  panel:  q  —  4.1295  a.u.,  a  =  49.3 
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the  fraction  of  the  unit  probability  associated  with  that  interval.  Typical  dis¬ 
tributions  for  q  =  0.0413  a.u.  (X  =  0.001),  q  =  1.8506  a.u.  (X  =  0.300),  and 
q  =  4.1295  a.u.  (X  =  1.000)  are  given  in  figure  8  and  these  distributions  are 
compared  graphically  with  the  normal  (Gaussian)  distribution,  having  the  same 
mean  and  variance  (that  is,  a2)  and  (6/2)  exp  (— \Az\/b)d(Az),  where  a  —  V2 b, 
which  is  the  double  exponential  distribution  used  by  Kendall  [7]  in  some  of  his 
theoretical  work. 
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